<html>
 <head>
  <meta charset="utf-8"/>
  <meta content="width=device-width, initial-scale=1, maximum-scale=1, user-scalable=no" name="viewport"/>
  <title>
   主题：求助，曲线间差异显著性分析  （简单一个问题为什么找不到解答）  | 数螺 | NAUT IDEA
  </title>
  <link href="http://cdn.bootcss.com/bootstrap/3.3.6/css/bootstrap-theme.min.css" rel="stylesheet"/>
  <link href="http://cdn.bootcss.com/bootstrap/3.3.6/css/bootstrap.min.css" rel="stylesheet"/>
  <style type="text/css">
   #xmain img {
                  max-width: 100%;
                  display: block;
                  margin-top: 10px;
                  margin-bottom: 10px;
                }

                #xmain p {
                    line-height:150%;
                    font-size: 16px;
                    margin-top: 20px;
                }

                #xmain h2 {
                    font-size: 24px;
                }

                #xmain h3 {
                    font-size: 20px;
                }

                #xmain h4 {
                    font-size: 18px;
                }


                .header {
	           background-color: #0099ff;
	           color: #ffffff;
	           margin-bottom: 20px;
	        }

	        .header p {
                  margin: 0px;
                  padding: 10px 0;
                  display: inline-block;  
                  vertical-align: middle;
                  font-size: 16px;
               }

               .header a {
                 color: white;
               }

              .header img {
                 height: 25px;
              }
  </style>
  <script src="http://cdn.bootcss.com/jquery/3.0.0/jquery.min.js">
  </script>
  <script src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML" type="text/javascript">
   MathJax.Hub.Config({elements: ["bbpress-forums"]});
  </script>
  <script src="http://nautstatic-10007657.file.myqcloud.com/static/css/readability.min.js" type="text/javascript">
  </script>
  <script type="text/javascript">
   $(document).ready(function() {
                 var loc = document.location;
                 var uri = {
                  spec: "http://cos.name/cn/topic/9912/",
                  host: "http://cos.name",
                  prePath: "http://cos.name",
                  scheme: "http",
                  pathBase: "http://cos.name/"
                 };
    
                 var documentClone = document.cloneNode(true);
                 var article = new Readability(uri, documentClone).parse();
     
                 document.getElementById("xmain").innerHTML = article.content;
                });
  </script>
  <!-- 1466438944: Accept with keywords: (title(0.333333333333):曲线,显著性,差异,主题,问题,论坛, topn(0.266666666667):曲线,帖子,主题,会员,显著性,题目,平均值,交叉,模型,例子,数据,系数,变量,效果,差异,问题,斜率,普通,用户名,随机,密码,笼统,结果,多多指教,显著,命令行,论坛,昆虫,连续变量,讨论区).-->
 </head>
 <body class="topic bbpress single single-topic postid-9912 single-author sidebar" onload="">
  <div class="header">
   <div class="container">
    <div class="row">
     <div class="col-xs-6 col-sm-6 text-left">
      <a href="/databee">
       <img src="http://nautidea-10007657.cos.myqcloud.com/logo_white.png"/>
      </a>
      <a href="/databee">
       <p>
        数螺
       </p>
      </a>
     </div>
     <div class="hidden-xs col-sm-6 text-right">
      <p>
       致力于数据科学的推广和知识传播
      </p>
     </div>
    </div>
   </div>
  </div>
  <div class="container text-center">
   <h1>
    主题：求助，曲线间差异显著性分析  （简单一个问题为什么找不到解答）
   </h1>
  </div>
  <div class="container" id="xmain">
   <div class="hfeed site" id="page">
    <header class="site-header" id="masthead" role="banner">
     <div id="cos-logo">
      <a href="http://cos.name/cn">
       <img src="http://cos.name/cn/wp-content/themes/COS-forest/images/headers/cos-logo.png"/>
      </a>
     </div>
     <div class="navbar" id="navbar">
      <nav class="navigation main-navigation" id="site-navigation" role="navigation">
       <h3 class="menu-toggle">
        菜单
       </h3>
       <div class="menu-%e8%8f%9c%e5%8d%951-container">
        <ul class="nav-menu" id="menu-%e8%8f%9c%e5%8d%951">
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-home menu-item-407772" id="menu-item-407772">
          <a href="http://cos.name/cn/">
           论坛首页
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-407773" id="menu-item-407773">
          <a href="http://cos.name/cn/forums/">
           讨论区
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-407774" id="menu-item-407774">
          <a href="http://cos.name/cn/wp-login.php?action=register">
           注册
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-407819" id="menu-item-407819">
          <a href="http://cos.name/">
           主站
          </a>
         </li>
        </ul>
       </div>
      </nav>
      <!-- #site-navigation -->
     </div>
     <!-- #navbar -->
    </header>
    <!-- #masthead -->
    <div class="site-main" id="main">
     <div class="content-area" id="primary">
      <div class="site-content" id="content" role="main">
       <article class="post-9912 topic type-topic status-publish hentry" id="post-9912">
        <header class="entry-header">
         <h1 class="entry-title">
          求助，曲线间差异显著性分析  （简单一个问题为什么找不到解答）
         </h1>
        </header>
        <!-- .entry-header -->
        <div class="entry-content">
         <div id="bbpress-forums">
          <div class="bbp-breadcrumb">
           <p>
            <a class="bbp-breadcrumb-home" href="http://cos.name/cn/">
             COS论坛 | 统计之都
            </a>
            <span class="bbp-breadcrumb-sep">
             ›
            </span>
            <a class="bbp-breadcrumb-root" href="http://cos.name/cn/forums/">
             讨论区
            </a>
            <span class="bbp-breadcrumb-sep">
             ›
            </span>
            <a class="bbp-breadcrumb-forum" href="http://cos.name/cn/forum/software/">
             软件应用
            </a>
            <span class="bbp-breadcrumb-sep">
             ›
            </span>
            <a class="bbp-breadcrumb-forum" href="http://cos.name/cn/forum/software/sas/">
             SAS
            </a>
            <span class="bbp-breadcrumb-sep">
             ›
            </span>
            <span class="bbp-breadcrumb-current">
             求助，曲线间差异显著性分析  （简单一个问题为什么找不到解答）
            </span>
           </p>
          </div>
          <div class="bbp-template-notice info">
           <p class="bbp-topic-description">
            该主题包含 9 条回复，3个帖子，最后由
            <a class="bbp-author-avatar" href="http://cos.name/cn/profile/32827/" rel="nofollow" title="查看jimi9876的档案">
             <img src="http://sdn.geekzu.org/avatar/f1698224cbfc3aea37c59355c49516f4?s=14&amp;d=monsterid&amp;r=g"/>
            </a>
            <a class="bbp-author-name" href="http://cos.name/cn/profile/32827/" rel="nofollow" title="查看jimi9876的档案">
             jimi9876
            </a>
            在
            <a href="http://cos.name/cn/topic/9912/#post-245074" title="回复：求助，曲线间差异显著性分析  （简单一个问题为什么找不到解答）">
             8 年, 3 月 之前
            </a>
            更新。
           </p>
          </div>
          <div class="bbp-pagination">
           <div class="bbp-pagination-count">
            查看 10 个帖子 - 1 到 10（总计 10 个）
           </div>
           <div class="bbp-pagination-links">
           </div>
          </div>
          <ul class="forums bbp-replies" id="topic-9912-replies">
           <li class="bbp-header">
            <div class="bbp-reply-author">
             作者
            </div>
            <!-- .bbp-reply-author -->
            <div class="bbp-reply-content">
             帖子
            </div>
            <!-- .bbp-reply-content -->
           </li>
           <!-- .bbp-header -->
           <li class="bbp-body">
            <div class="bbp-reply-header" id="post-9912">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2008年3月10日 下午3:58
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/9912/#post-9912">
               1 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-9912 -->
            <div class="odd bbp-parent-forum-998 bbp-parent-topic-9912 bbp-reply-position-1 user-id-32827 topic-author post-9912 topic type-topic status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/32827/" rel="nofollow" title="查看jimi9876的档案">
               <img src="http://sdn.geekzu.org/avatar/f1698224cbfc3aea37c59355c49516f4?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/32827/" rel="nofollow" title="查看jimi9876的档案">
               jimi9876
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <p>
               这是一个经常碰到的问题，我查了网上以及参考书似乎都没有相似的例子或权威的回答，还请您不吝赐教，将十分感谢！
              </p>
              <p>
               本数据是一种昆虫经三种处理后测定了3天的产卵情况（每个处理3个重复），如果用每个处理每一天产卵数的均值相对于天数作图，将是三条下降的曲线，请问如何检测三条曲线（三个处理三天内产卵下降趋势）差异显著性（并作成对比较），原始数据如下
              </p>
              <p>
               Treat    Day    Response
               <br/>
               T1    1    15.58
               <br/>
               T1    1    14.69
               <br/>
               T1    1    14.01
               <br/>
               T1    2    13.01
               <br/>
               T1    2    12.12
               <br/>
               T1    2    11.43
               <br/>
               T1    3    10.8
               <br/>
               T1    3    10.2
               <br/>
               T1    3    9.825
               <br/>
               T2    1    15.313
               <br/>
               T2    1    13.22
               <br/>
               T2    1    12.04
               <br/>
               T2    2    10.833
               <br/>
               T2    2    9.753
               <br/>
               T2    2    8.973
               <br/>
               T2    3    8.277
               <br/>
               T2    3    7.167
               <br/>
               T2    3    6.669
               <br/>
               T3    1    11.313
               <br/>
               T3    1    9.22
               <br/>
               T3    1    8.04
               <br/>
               T3    2    6.833
               <br/>
               T3    2    5.753
               <br/>
               T3    2    4.973
               <br/>
               T3    3    4.277
               <br/>
               T3    3    3.167
               <br/>
               T3    3    2.669
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-244519">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2008年3月11日 上午4:16
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/9912/#post-244519">
               2 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-244519 -->
            <div class="even bbp-parent-forum-998 bbp-parent-topic-9912 bbp-reply-position-2 user-id-22116 post-244519 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/22116/" rel="nofollow" title="查看hexm26的档案">
               <img src="http://sdn.geekzu.org/avatar/2d9c052731bdfc1c7f4b06e84fe8d2f4?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/22116/" rel="nofollow" title="查看hexm26的档案">
               hexm26
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <p>
               信息不全！你的第一行的15.58和第四行的13.01来自同一个昆虫吗？
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-244532">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2008年3月11日 上午5:02
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/9912/#post-244532">
               3 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-244532 -->
            <div class="odd bbp-parent-forum-998 bbp-parent-topic-9912 bbp-reply-position-3 user-id-13656 post-244532 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/13656/" rel="nofollow" title="查看王笑权的档案">
               <img src="http://sdn.geekzu.org/avatar/d01c79fb81e8c445185b8d27ac03577f?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/13656/" rel="nofollow" title="查看王笑权的档案">
               王笑权
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <p>
               不喝完8两老白干儿解决不了这个问题。 干杯！
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-244562">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2008年3月11日 上午11:27
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/9912/#post-244562">
               4 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-244562 -->
            <div class="even bbp-parent-forum-998 bbp-parent-topic-9912 bbp-reply-position-4 user-id-32827 topic-author post-244562 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/32827/" rel="nofollow" title="查看jimi9876的档案">
               <img src="http://sdn.geekzu.org/avatar/f1698224cbfc3aea37c59355c49516f4?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/32827/" rel="nofollow" title="查看jimi9876的档案">
               jimi9876
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <blockquote class="d4pbbc-quote">
               <p>
                [b]引用第1楼[i]hexm26[/i]于[i]2008-03-11 11:16[/i]发表的“”[/b]:
                <br/>
                信息不全！你的第一行的15.58和第四行的13.01来自同一个昆虫吗？
               </p>
              </blockquote>
              <p>
               这回前了，Insect一列是昆虫编号,谢谢，还请多多指教！
              </p>
              <p>
               Treat    Day    Insect    Response
               <br/>
               T1    1    1    15.58
               <br/>
               T1    1    2    14.69
               <br/>
               T1    1    3    14.01
               <br/>
               T1    2    1    13.01
               <br/>
               T1    2    2    12.12
               <br/>
               T1    2    3    11.43
               <br/>
               T1    3    1    10.8
               <br/>
               T1    3    2    10.2
               <br/>
               T1    3    3    9.825
               <br/>
               T2    1    4    15.313
               <br/>
               T2    1    5    13.22
               <br/>
               T2    1    6    12.04
               <br/>
               T2    2    4    10.833
               <br/>
               T2    2    5    9.753
               <br/>
               T2    2    6    8.973
               <br/>
               T2    3    4    8.277
               <br/>
               T2    3    5    7.167
               <br/>
               T2    3    6    6.669
               <br/>
               T3    1    7    11.313
               <br/>
               T3    1    8    9.22
               <br/>
               T3    1    9    8.04
               <br/>
               T3    2    7    6.833
               <br/>
               T3    2    8    5.753
               <br/>
               T3    2    9    4.973
               <br/>
               T3    3    7    4.277
               <br/>
               T3    3    8    3.167
               <br/>
               T3    3    9    2.669
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-244585">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2008年3月12日 上午2:12
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/9912/#post-244585">
               5 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-244585 -->
            <div class="odd bbp-parent-forum-998 bbp-parent-topic-9912 bbp-reply-position-5 user-id-22116 post-244585 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/22116/" rel="nofollow" title="查看hexm26的档案">
               <img src="http://sdn.geekzu.org/avatar/2d9c052731bdfc1c7f4b06e84fe8d2f4?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/22116/" rel="nofollow" title="查看hexm26的档案">
               hexm26
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <p>
               用随机系数模型（Random Coefficients Model）吧：
               <br/>
               data insect;
               <br/>
               input Treat $ Day Insect Response;
               <br/>
               datalines;
               <br/>
               T1    1    1    15.58
               <br/>
               T1    1    2    14.69
               <br/>
               T1    1    3    14.01
               <br/>
               T1    2    1    13.01
               <br/>
               T1    2    2    12.12
               <br/>
               T1    2    3    11.43
               <br/>
               T1    3    1    10.8
               <br/>
               T1    3    2    10.2
               <br/>
               T1    3    3    9.825
               <br/>
               T2    1    4    15.313
               <br/>
               T2    1    5    13.22
               <br/>
               T2    1    6    12.04
               <br/>
               T2    2    4    10.833
               <br/>
               T2    2    5    9.753
               <br/>
               T2    2    6    8.973
               <br/>
               T2    3    4    8.277
               <br/>
               T2    3    5    7.167
               <br/>
               T2    3    6    6.669
               <br/>
               T3    1    7    11.313
               <br/>
               T3    1    8    9.22
               <br/>
               T3    1    9    8.04
               <br/>
               T3    2    7    6.833
               <br/>
               T3    2    8    5.753
               <br/>
               T3    2    9    4.973
               <br/>
               T3    3    7    4.277
               <br/>
               T3    3    8    3.167
               <br/>
               T3    3    9    2.669
               <br/>
               ;;;;
              </p>
              <p>
               proc mixed data=insect;
               <br/>
               class treat insect;
               <br/>
               model Response = treat day treat*day;
               <br/>
               random intercept day / type=un sub=insect;
               <br/>
               lsmeans treat / adjust = TUKEY ;
               <br/>
               run;
              </p>
              <p>
               结论：三条曲线之间差异显著（P=0.0200），其中T3显著差异于T1（P=0.0001）和T2（P=0.0023），T1和T2之间也存在差异但显著性在10％程度内（P=0.0658）。
              </p>
              <p>
               上面的分析并不唯一，用重复测量（Repeated Measures）或其他的模型（type=ar(1)）等也可以做，但结论应该差不多。另外，我用的的Tukey调整多元比较，具体情况楼主自己根据试验目的调整吧。
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-244799">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2008年3月13日 上午12:00
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/9912/#post-244799">
               6 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-244799 -->
            <div class="even bbp-parent-forum-998 bbp-parent-topic-9912 bbp-reply-position-6 user-id-32827 topic-author post-244799 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/32827/" rel="nofollow" title="查看jimi9876的档案">
               <img src="http://sdn.geekzu.org/avatar/f1698224cbfc3aea37c59355c49516f4?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/32827/" rel="nofollow" title="查看jimi9876的档案">
               jimi9876
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <blockquote class="d4pbbc-quote">
               <p>
                [b]引用第4楼[i]hexm26[/i]于[i]2008-03-12 09:12[/i]发表的“”[/b]:
                <br/>
                用随机系数模型（Random Coefficients Model）吧：
                <br/>
                data insect;
                <br/>
                input Treat $ Day Insect Response;
                <br/>
                datalines;
                <br/>
                T1    1    1    15.58
                <br/>
                …….
               </p>
              </blockquote>
              <p>
               非常感谢！虽然我不懂sas，但用mixed 或repeat模型应该是对的，但用lsmeans treat / adjust = TUKEY做成对比较似乎有问题，比如我对数据做了调整，让T1与T3相同，即随时间下降，而T2则随时间上升：
               <br/>
               Treat    Day    Insect    Response
               <br/>
               T1    1    1    15.58
               <br/>
               T1    1    2    14.69
               <br/>
               T1    1    3    14.01
               <br/>
               T1    2    1    13.01
               <br/>
               T1    2    2    12.12
               <br/>
               T1    2    3    11.43
               <br/>
               T1    3    1    10.8
               <br/>
               T1    3    2    10.2
               <br/>
               T1    3    3    9.825
               <br/>
               T2    1    4    9.825
               <br/>
               T2    1    5    10.2
               <br/>
               T2    1    6    10.8
               <br/>
               T2    2    4    11.43
               <br/>
               T2    2    5    12.12
               <br/>
               T2    2    6    13.01
               <br/>
               T2    3    4    14.01
               <br/>
               T2    3    5    14.69
               <br/>
               T2    3    6    15.58
               <br/>
               T3    1    7    15.58
               <br/>
               T3    1    8    14.69
               <br/>
               T3    1    9    14.01
               <br/>
               T3    2    7    13.01
               <br/>
               T3    2    8    12.12
               <br/>
               T3    2    9    11.43
               <br/>
               T3    3    7    10.8
               <br/>
               T3    3    8    10.2
               <br/>
               T3    3    9    9.825
              </p>
              <p>
               则结果为
               <br/>
               Null Model Likelihood Ratio Test
              </p>
              <p>
               DF    Chi-Square      Pr &gt; ChiSq
              </p>
              <p>
               2         22.47          &lt;.0001
              </p>
              <p>
               Type 3 Tests of Fixed Effects
              </p>
              <p>
               Num     Den
               <br/>
               Effect         DF      DF    F Value    Pr &gt; F
              </p>
              <p>
               Treat           2       9      94.84    &lt;.0001
               <br/>
               Day             1       6     148.35    &lt;.0001
               <br/>
               Day*Treat       2       9     593.40    &lt;.0001
              </p>
              <p>
               Least Squares Means
              </p>
              <p>
               Standard
               <br/>
               Effect    Treat    Estimate       Error      DF    t Value    Pr &gt; |t|
              </p>
              <p>
               Treat     T1        12.4072      0.4100       9      30.26      &lt;.0001
               <br/>
               Treat     T2        12.4072      0.4100       9      30.26      &lt;.0001
               <br/>
               Treat     T3        12.4072      0.4100       9      30.26      &lt;.0001
              </p>
              <p>
               Effect   Treat   _Treat   Estimate      Error     DF   t Value   Pr &gt; |t|   Adjustment      Adj P
              </p>
              <p>
               Treat    T1      T2       6.52E-15     0.5798      9      0.00     1.0000   Tukey          1.0000
               <br/>
               Treat    T1      T3       -947E-17     0.5798      9     -0.00     1.0000   Tukey          1.0000
               <br/>
               Treat    T2      T3        -16E-15     0.5798      9     -0.00     1.0000   Tukey          1.0000
              </p>
              <p>
               此成对比较显然与上述结果矛盾，还请您多多指教，十分感谢！
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-244809">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2008年3月13日 上午2:13
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/9912/#post-244809">
               7 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-244809 -->
            <div class="odd bbp-parent-forum-998 bbp-parent-topic-9912 bbp-reply-position-7 user-id-22116 post-244809 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/22116/" rel="nofollow" title="查看hexm26的档案">
               <img src="http://sdn.geekzu.org/avatar/2d9c052731bdfc1c7f4b06e84fe8d2f4?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/22116/" rel="nofollow" title="查看hexm26的档案">
               hexm26
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <p>
               原因很简单，你的问题开始说的是3条曲线都是下降的，所以我给你的模型中Day是个连续变量，用lsmeans treat / adjust = TUKEY其实比较的是一个平均值。你刻意更改的数据明显是3条曲线的平均值完全想同。
               <br/>
               既然你改了题目，分析也要相应的改动，设置一个新的变量Dayc=Day，但Dayc是离散变量，将它放在Class和Model命令行里。然后用Contrast或Estimate分析每天的曲线差异。记得SAS的Help文件里面有个类似的例子。我现在手边没有SAS，你自己改改把结果贴上来，我们一起来分析。如果我的估计没错的话，3天T2和T1（T3）的差别有，但都不显著。
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-244987">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2008年3月14日 上午11:34
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/9912/#post-244987">
               8 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-244987 -->
            <div class="even bbp-parent-forum-998 bbp-parent-topic-9912 bbp-reply-position-8 user-id-32827 topic-author post-244987 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/32827/" rel="nofollow" title="查看jimi9876的档案">
               <img src="http://sdn.geekzu.org/avatar/f1698224cbfc3aea37c59355c49516f4?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/32827/" rel="nofollow" title="查看jimi9876的档案">
               jimi9876
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <blockquote class="d4pbbc-quote">
               <p>
                [b]引用第6楼[i]hexm26[/i]于[i]2008-03-13 09:13[/i]发表的“”[/b]:
                <br/>
                原因很简单，你的问题开始说的是3条曲线都是下降的，所以我给你的模型中Day是个连续变量，用lsmeans treat / adjust = TUKEY其实比较的是一个平均值。你刻意更改的数据明显是3条曲线的平均值完全想同。
                <br/>
                既然你改了题目，分析也要相应的改动，设置一个新的变量Dayc=Day，但Dayc是离散变量，将它放在Class和Model命令行里。然后用Contrast或Estimate分析每天的曲线差异。记得SAS的Help文件里面有个类似的例子。我现在手边没有SAS，你自己改改把结果贴上来，我们一起来分析。如果我的估计没错的话，3天T2和T1（T3）的差别有，但都不显著。
               </p>
              </blockquote>
              <p>
               非常感谢，实在非常感谢！看了不少书，确实是用contrast两两比较，编了个程序，还请指正
              </p>
              <p>
               data insect;
               <br/>
               input Treat $ Day Insect Response;
               <br/>
               datalines;
               <br/>
               T1   1   1   15.58
               <br/>
               T1   1   2   14.69
               <br/>
               T1   1   3   14.01
               <br/>
               T1   2   1   13.01
               <br/>
               T1   2   2   12.12
               <br/>
               T1   2   3   11.43
               <br/>
               T1   3   1   10.8
               <br/>
               T1   3   2   10.2
               <br/>
               T1   3   3   9.825
               <br/>
               T2   1   4   15.313
               <br/>
               T2   1   5   13.22
               <br/>
               T2   1   6   12.04
               <br/>
               T2   2   4   10.833
               <br/>
               T2   2   5   9.753
               <br/>
               T2   2   6   8.973
               <br/>
               T2   3   4   8.277
               <br/>
               T2   3   5   7.167
               <br/>
               T2   3   6   6.669
               <br/>
               T3   1   7   11.313
               <br/>
               T3   1   8   9.22
               <br/>
               T3   1   9   8.04
               <br/>
               T3   2   7   6.833
               <br/>
               T3   2   8   5.753
               <br/>
               T3   2   9   4.973
               <br/>
               T3   3   7   4.277
               <br/>
               T3   3   8   3.167
               <br/>
               T3   3   9   2.669
               <br/>
               ;;;;
              </p>
              <p>
               proc mixed data=insect;
               <br/>
               class treat day insect;
               <br/>
               model Response = treat day treat*day;
               <br/>
               repeated day /  sub=insect type=un r rcorr;
               <br/>
               contrast 'treat T1  vs T2 intercept' treat 1 -1 0;
               <br/>
               contrast 'treat T1  vs T2 slope' day 1  -1 0;
              </p>
              <p>
               run;
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-245047">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2008年3月15日 上午4:06
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/9912/#post-245047">
               9 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-245047 -->
            <div class="odd bbp-parent-forum-998 bbp-parent-topic-9912 bbp-reply-position-9 user-id-22116 post-245047 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/22116/" rel="nofollow" title="查看hexm26的档案">
               <img src="http://sdn.geekzu.org/avatar/2d9c052731bdfc1c7f4b06e84fe8d2f4?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/22116/" rel="nofollow" title="查看hexm26的档案">
               hexm26
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <p>
               嗯,不错,你用Repeated Measures Model来做也可以.另外,treat*day这个交叉如果是有显著作用的话,你在分析treat的主成分效果时应该把交叉项也带上,也就是说每天的曲线效果是不同的.在曲线有明显不同斜率的情况下,笼统的说treat的总效果是没有意义的.具体如何实现,你可以通过Contrast来做,比如说treat 1 -1 0 treat*day 1 0 0 -1 0 0 0 0 0.不一定准确,大概是这个样子的,你自己去试吧.
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-245074">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2008年3月15日 上午10:30
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/9912/#post-245074">
               10 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-245074 -->
            <div class="even bbp-parent-forum-998 bbp-parent-topic-9912 bbp-reply-position-10 user-id-32827 topic-author post-245074 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/32827/" rel="nofollow" title="查看jimi9876的档案">
               <img src="http://sdn.geekzu.org/avatar/f1698224cbfc3aea37c59355c49516f4?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/32827/" rel="nofollow" title="查看jimi9876的档案">
               jimi9876
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <blockquote class="d4pbbc-quote">
               <p>
                [b]引用第8楼[i]hexm26[/i]于[i]2008-03-15 11:06[/i]发表的“”[/b]:
                <br/>
                嗯,不错,你用Repeated Measures Model来做也可以.另外,treat*day这个交叉如果是有显著作用的话,你在分析treat的主成分效果时应该把交叉项也带上,也就是说每天的曲线效果是不同的.在曲线有明显不同斜率的情况下,笼统的说treat的总效果是没有意义的.具体如何实现,你可以通过Contrast来做,比如说treat 1 -1 0 treat*day 1 0 0 -1 0 0 0 0 0.不一定准确,大概是这个样子的,你自己去试吧.
               </p>
              </blockquote>
              <p>
               真是学无止境，非常感谢您的决定性指导！
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
           </li>
           <!-- .bbp-body -->
           <li class="bbp-footer">
            <div class="bbp-reply-author">
             作者
            </div>
            <div class="bbp-reply-content">
             帖子
            </div>
            <!-- .bbp-reply-content -->
           </li>
           <!-- .bbp-footer -->
          </ul>
          <!-- #topic-9912-replies -->
          <div class="bbp-pagination">
           <div class="bbp-pagination-count">
            查看 10 个帖子 - 1 到 10（总计 10 个）
           </div>
           <div class="bbp-pagination-links">
           </div>
          </div>
          <div class="bbp-no-reply" id="no-reply-9912">
           <div class="bbp-template-notice">
            <p>
             您必须先登录才能回复该主题。
            </p>
           </div>
          </div>
         </div>
        </div>
        <!-- .entry-content -->
        <footer class="entry-meta">
        </footer>
        <!-- .entry-meta -->
       </article>
       <!-- #post -->
       <div class="comments-area" id="comments">
       </div>
       <!-- #comments -->
      </div>
      <!-- #content -->
     </div>
     <!-- #primary -->
     <div class="sidebar-container" id="tertiary" role="complementary">
      <div class="sidebar-inner">
       <div class="widget-area">
        <aside class="widget bbp_widget_login" id="bbp_login_widget-2">
         <h3 class="widget-title">
          登录
         </h3>
         <form action="http://cos.name/cn/wp-login.php" class="bbp-login-form" method="post">
          <fieldset>
           <legend>
            登录
           </legend>
           <div class="bbp-username">
            <label for="user_login">
             用户名:
            </label>
           </div>
           <div class="bbp-password">
            <label for="user_pass">
             密码:
            </label>
           </div>
           <div class="bbp-remember-me">
            <label for="rememberme">
             记住用户名
            </label>
           </div>
           <div class="bbp-submit-wrapper">
            <button class="button submit user-submit" id="user-submit" name="user-submit" tabindex="104" type="submit">
             登录
            </button>
           </div>
           <div class="bbp-login-links">
            <a class="bbp-register-link" href="http://cos.name/cn/wp-login.php?action=register" title="注册">
             注册
            </a>
            <a class="bbp-lostpass-link" href="http://cos.name/cn/wp-login.php?action=lostpassword" title="忘记密码">
             忘记密码
            </a>
           </div>
          </fieldset>
         </form>
        </aside>
        <aside class="widget widget_text" id="text-7">
         <h3 class="widget-title">
          搜索
         </h3>
         <div class="textwidget">
          <form action="http://www.google.com/search" id="bbp-search-form" method="get" onsubmit="Gsitesearch(this)" role="search">
           <div>
           </div>
          </form>
          <form id="bbp-search-form-baidu" onsubmit="g(this)" role="search">
           <div>
           </div>
          </form>
         </div>
        </aside>
        <aside class="widget widget_text" id="text-2">
         <h3 class="widget-title">
          新鲜事
         </h3>
         <div class="textwidget">
          <ul>
           <li>
            <a href="http://cos.name/cn/topics/">
             最新帖子
            </a>
           </li>
           <li>
            <a href="http://cos.name/cn/view/popular/">
             最热门主题
            </a>
           </li>
           <li>
            <a href="http://cos.name/cn/view/no-replies/">
             消灭零回复
            </a>
           </li>
          </ul>
         </div>
        </aside>
        <aside class="widget widget_text" id="text-3">
         <h3 class="widget-title">
          RSS订阅
         </h3>
         <div class="textwidget">
          <ul>
           <li>
            <img src="http://cos.name/wp-includes/images/rss.png"/>
            <a href="http://cos.name/cn/topics/feed/">
             所有主题
            </a>
           </li>
           <li>
            <img src="http://cos.name/wp-includes/images/rss.png"/>
            <a href="http://cos.name/cn/forums/feed/">
             所有帖子
            </a>
           </li>
          </ul>
         </div>
        </aside>
       </div>
       <!-- .widget-area -->
      </div>
      <!-- .sidebar-inner -->
     </div>
     <!-- #tertiary -->
    </div>
    <!-- #main -->
    <footer class="site-footer" id="colophon" role="contentinfo">
     <div class="site-info">
      版权所有 © 2014 统计之都 | 由
      <a href="http://wordpress.org/">
       WordPress
      </a>
      构建 | 主题修改自
      <a href="http://wordpress.org/themes/twentythirteen">
       Twenty Thirteen
      </a>
     </div>
     <!-- .site-info -->
    </footer>
    <!-- #colophon -->
   </div>
   <!-- #page -->
  </div>
 </body>
</html>